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ABSTRACT 

Context. Since the early 1990s Gamma Ray Bursts have been accepted to be of extra-galactic origin due to the isotropic distribution 
observed by BATSE and the redshifts observed via absorption line spectroscopy. Nevertheless, upon further examination at least one 
case turned out to be of galactic origin. This particular event presented a Fast Rise, Exponential Decay (FRED) structure which leads 
us to believe that other FRED sources might also be Galactic. 
Aims. This study was set out to estimate the most probable degree of contamination by galactic sources that certain samples of FREDs 
have. 



in 



Methods. In order to quantify the degree of anisotropy the average dipolar and quadripolar moments of each sample of GRBs with 
respect to the galactic plane were calculated. This was then compared to the probability distribution of simulated samples comprised 
of a combination of isotropically generated sources and galactic sources. 

Results. We observe that the dipolar and quadripolar moments of the selected subsamples of FREDs are found more than two standard 
deviations outside those of random isotropically generated samples.The most probable degree of contamination by galactic sources 
for the FRED GRBs of the Swift catalog detected until February 201 1 that do not have a known redshift is about 21 out of 77 sources 
which representes roughly equal to 27%. Furthermore we observe, that by removing from this sample those bursts that have any type 
of indirect redshift indicator and multiple peaks the most probable contamination increases up to 34% (17 out of 49 sources). 
Conclusions. It is probable that a high degree of contamination by galactic sources occurs among the single peak FREDs observed 
by Swift. So we encourage further studies on these type of events in order to determine the nature of what could be an exotic type of 
\ galactic source. 

^\ Key words, (stars:) Gamma-ray burst: general - stars: magnetars - methods: statistical 

ON 
VO, 

^1 ■ 1. Introduction emission is emitted, where one or multiple peaks with a vari- 
ety of shapes are observed. However, some of them present a 

I Gamma Ray Bursts (GRBs) are the most energetic explosions Fast Rise and Exponential Decay (FRED hereafter) behavior and 

(N ■ known in the Universe second only to the Big Bang. Discovered these have been correlated with other properties of the bursts 

r-H ■ in the 1960s, they were widely believed to originate in the Milky (|Bhat et al.ll 19941) . suggesting that they may be of a different na- 

^ Way because of their relatively high flux of photons which needs ture to that f other GRBs. 
an unprecedented emission mechanism to account for such a 

high energy output. It was not until 1997 when the first measure- Thgre has been at leafjt one reported GRB mat upon fur . 

H . ment of redshift was performed on a GRB afterglow that it was ther examinati on has resulted to be phenomena from wit hin the 

asserted without doubt the cosmological nature of these objects Milky Way dCastro-Tirado et al.ll2008l IStefanescu et al.l l2008). 

fl Metzger et al. || 1997| ). This source displayed a Fast Rise, Exponential Decay (FRED) 

GRB afterglows fade within a few hours and as a conse- structure which leads to believe that there could be others like it. 
quence the redshift of most GRBs are unknown. In the past 

several studies have been carried out to indirectly determine the In mis wor k, we w m try to estimate the most probable de- 

galactic or extragalactic nature of the bursts by analyzing their gree of contamination by galactic sources that certain samples 

spatial distribution in the sky (|Balazs et al.||l998t |Mazets et aj| D f FREDs have. We have organized the paper as follows: ir 

119811 |Meegan et al.|ll992j), and historically it ser ved as a strong Section 2 we establish the selection criteria of the studied sam 

argument against the galactic origin of GRBs (|Paciesas et al.| p i es . Section 3 describes the methodology used for quantify- 

|1999|). This technique has also been used to suggest a more local ing the anisotropy and determining the probability of observing 

nature of long-lag bursts by showing that they may be relate d to these values for both extragalactic and galactic sources while 

the Super-Galactic structure ([Foley et al.||2008t |Norris||2002|). taking in to account the exposure of Swift. We discuss the results 

The observed light curve of each GRB varies from burst to from our analysis in Section 4, and give main conclusions in 

burst, specially during in the prompt phase when the gamma ray Section 5. 
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2. Sample selection 

In order to have a homogeneous distribution, only Swi/f-detected 
GRBs were taken into account. From the catalog of 596 GRBs 
detected by Swift before March 201 1, 111 of the GRBs were se- 
lected for having a FRED structure reported in a GCN. Using the 
information available in peer reviewed paper^] and other GCN 
circulars, related to the 111 FRED GRBs, the following subsam- 
ples were selected: 

- Sample 1: All 111 FREDs detected by Swift up until 
February 2011. 

- Sample 2: 77 FREDs without any measured redshift. 

- Sample 3: 71 FREDs without stated "high redshift criteria". 

- Sample 4: 59 FREDs without any type of indirect redshift 
indication. 

- Sample 5: 49 FREDs without any redshift indications or 
multiple peaks. 




Sample 1 • 


Sample 3 • 


Sample 5 • 


Sample 2 • 


Sample 4 





Fig. 1. Sources in each one of the samples analized. Note that to 
avoid redundancy only the sources not present in the subsequent 
samples were included. 

It is important to note that only Sample 5 included solely 
those bursts that consisted of one pure FRED peak. 

3. Anisotropy quantification 

It has been proven ( Hartmann & Epstein 1989) that the mean 
dipolar and quadripolar moments of the galactic coordinates 
(cosb and sin 2 b respectively, where b is the galactic latitude) 
are good t ools to quantify the isotropy with respect to the galac- 
tic plane ([Castro Tiradolll994h . The degree of isotropy of each 
sample was calculated using the coordinates available from the 
GCNs for each burst and the results are shown in Table Q] 

3.1. Exposure map 

Due to the nature of its instruments, orbit and mission, Swift's 
pointing towards the sky is not homogeneous. It is of particu- 
lar relevance to note that there has been less integrated expo- 
sure time towards the Milky Way's disk than toward the galac- 
tic poles. This fact would represent a bias for the nature of the 

1 Only two peer reviewed papers affected the sample selection 
dClemens et al.ll201 lhlPeriev et all2009h 

2 For a list of which specific bursts where selected see the arXiv ver- 
sion of this paper: arXiv: 1210.3699 [astro-ph.HE] 



Sample 


< cosb > 


< sirrb > 


#1 


0.7883 


0.3397 


#2 


0.8221 


0.2860 


#3 


0.8184 


0.2909 


#4 


0.8344 


0.2673 


#5 


0.8397 


0.2622 



Table 1. The dipolar and quadripolar moments of the samples 
as a quantitative measurement of the degree of isotropy in the 
samples. 



study carried out for this publication if left unaccounted, so a 
map was created by integrating the exposure mask function for 
the BAT instrument, multiplied by the exposure times of all ob- 
servations carried out between April 16, 2005 and February, 1st 
2011, taking into account the pointing and rotation of the BAT 
instrument 




Fig. 2. Swift exposure map in Galactic Coordinates derived for 
this study. Colors represent the exposure time (in seconds). 



3.2. Montecarlo simulations 

Montecarlo simulations were carried out to determine the 
Probability Mass Function (PMF) of the average dipolar and 
quadripolar moments of random GRB distributions. Therefore 
random coordinates were generated taking care that they would 
have a homogeneous distribution on an spherical surface. 

These random points were then used to determine the PMF 
of the dipolar and quadripolar moments of random sources in the 
sky in order to determine how much the observed samples' val- 
ues deviated from those of a completely isotropically generated 
sample. This was done generating an equal number of random 
points to that of each sample, recording the value of the mean 
dipolar and quadripolar moment and iterating a statistically sig- 
nificant number of times(10 6 - 10 9 iterations), the histogram of 
the recorded values was then used to determine the values for 
standard deviations (<x, 2<r, 3cr). 

3.3. Metropolis-Hastings algorithm 

To account for the anisotropy of Swiff s exposure of the night 
sky, it was necessary to factor in the probability that a particular 
random source was detected by Swift. This was done by using 
the Metropolis-Hastings algorithm which can be summarized in 
three steps: 

1. Create a random source (set of coordinates). 

3 The method used to derive the exposure map is the same as the one 
detailed in Veres et ah, 2010 
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Sample 


Dipolar 


Quadripolar 


#1 


78.24 


13.26 


#2 


97.91 


0.59 


#3 


96.72 


1.03 


#4 


98.52 


0.44 


#5 


97.83 


0.39 



Table 2. Percentile of the values observed for the dipolar and 
quadripolar moments of each sample. 



2. Create a random number with a range equal to the values of 
the map or mask. 

3. Compare the value of the map at those coordinates to the 
random number. 

4. If the value of the point on the map exceeds the random num- 
ber, then the random source is included in the sample for fur- 
ther analysis. Otherwise a new random source is created and 
the process is repeated until the correct amount of sources 
are obtained. 

This will effectively generate random sources that are more 
likely to appear where the exposure is higher. This method was 
tested by generating a statistically significant number of random 
sources and checking that the resulting image was one propor- 
tional to the weighting mask well within normal statistical fluc- 
tuations. 

3.4. Contamination by galactic sources 

Considering that i) the density of matter of the Milky Way is 
roughly correlated with the amount of interstellar dust, and by 
consequence so is the amount of stellar sources as well, and 
ii) the transparency of gam ma-rays to interstella r dust, we used 
maps of dust IR emission ( Schle gel et al.l [T998) as a weighting 
mask for the Metropolis-Hastings algorithm to generate random 
galactic sources. 

The isotropically generated samples were contaminated by 
increasing the number of galactically generated random sources 
(N) to observe how this affected the PMF of their dipolar and 
quadripolar moment. This was done considering all possible 
combinations for the number of GRBs in the different samples 
and taking into account the Swift exposure map for each source 
generated. 



4. Results 

The Montecarlo simulations of the isotropically generated ran- 
dom samples (weighted by the Swift exposure map) showed that 
the dipolar and quadripolar moments from the real samples were 
consistently deviated from the average as is shown in Figure [3] 

In table|2]the percentile of the population in which each sam- 
ple is located is shown. Considering that by definition, one stan- 
dard deviation will be between the 15.87th and the 84.13th per- 
centiles, two standard deviations between 2.28th and 97.72th and 
three deviations between 0.13th and 99.87th, we observe that 
with the exception of the first sample, all samples have dipolar 
and quadripolar moments located outside two standard devia- 
tions. 

The probability distribution of samples that contained both 
isotropically and galactically generated sources allowed us to 
compare how the contamination by galactic sources affected 
the likelihood of obtaining certain momentum values, for ex- 
ample see figure [4] This technique is similar to the one used 




0.7 0.75 0.8 0.85 

Dipolar Moment of each Sample <cos b> 



0.9 




0.25 0.3 0.35 0.4 0.45 0.5 

Quadripolar Moment of each sample <sen 2 1» 



0.55 



Fig. 3. Dipolar (top) and quadripolar (bottom) moments PMFs 
for samples of isotropically generated sources weighted by 
Swift's exposure map for each sample size, and the observed val- 
ues for each sample (vertical lines). 



Sample #: 


1 


2 


3 


4 


5 


Dipolar: 


13 


20 


17 


18 


16 


Dipolar %: 


11.71% 


25.97% 


23.94% 


30.51% 


32.65% 


Quadripolar: 


12 


22 


20 


21 


18 


Quad. %: 


10.81% 


28.57% 


28.17% 


35.59% 


36.73% 



Table 3. Amount of galactic sources that yield the higher prob- 
ability to obtain the observed dipolar and quadripolar moments 
and the percentage of the sample size they imply. 



in the past to study the degree of contamination by galactic 
repeater gamma-ray s ources present in two GRB catalogues 
dGorosabel et al.|[l998l) . By observing the probability of the ob- 
served values in each one of the curves that resulted from the 
simulations, we determined the relative probability that each one 
of those combinations of isotropically and galactically generated 
sources would yield the observed momentums. In figure [5] the 
probability as a function of the amount of galactic sources intro- 
duced in each sample is shown. 

5. Conclusions 

It is observed that, with the exception of the first sample, all ob- 
served samples show dipolar and quadripolar moments outside 
two standard deviations from the mean of an isotropically gener- 
ated distribution. Although this result is not conclusive, there is a 
high probability that the samples are not of a purely extragalactic 
nature. 

The probability of obtaining the dipolar and quadripolar mo- 
ments that are measured in the samples is much higher when 
including a significant amount of galactic sources than it is for 
only isotropically generated sources. 

As shown in table[3] if we consider the amount of contamina- 
tion that yields the highest probability of obtaining the observed 
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Fig. 4. Example of dipolar (top) and quadripolar (bottom) mo- 
ment probability distributions for samples of 49 randomly gen- 
erated sources with an increasing number (N) of those sources 
being of galactic origin, and the observed value for Sample #5 
(dashed vertical line). Each line has a different amount of galac- 
tic sources, starting with N=0 (red), and only every fifth line was 
colored for an easier reading. 

values, the amount of galactic sources that are probably contam- 
inating the Swift GRB catalog is between 16 and 22. This value 
represents approximately 3% of the catalog used for this study. 

Sample 5 has been narrowed down so that it is likely that 
1 out of every 3 is in fact a Galactic source, so it is of great 
interest to further study these sources to determine if there are 
other indications that they are not GRBs. 

The high Galactic extinction persuades optical ground based 
spectroscopy of most low galactic latitude GRBs. With the 
present work we show that a large fraction of those abandoned 
follow-ups could reveal a missing population of galactic events. 
We encourage ground observers to follow up those events, since 
it might open to the discovery of unknown high-energy phenom- 
ena present in our Galaxy. 
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Fig. 5. Relative probability of obtaining the dipolar (top) and 
quadripolar (bottom) moments ±0.002 measured for our sam- 
ples. 
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This appendix specifies which bursts are included in each one of the five samples used in this study. To facilitate its reading each 
column was named in reference to the discriminating parameter for each subsample. That is: 

- Sample 1: All 1 1 1 FREDs detected by Swift before March 201 1. 

- Sample 2: 77 FREDs without any measured redshift different from zero 

- Sample 3: 71 FREDs without any claimed high redshift criteria 

- Sample 4: 59 FREDs without any type of indirect redshift indication 

- Sample 5: 49 FREDs without any redshift indications or multiple peaks 

With the exception of the first column and where noted otherwise, each number refers to the GCN circular. 
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GRB# 


F.R.E.D. 


Redshift 


"High-z" 


Any z 


Peaks 


071028B 


7019 


none 


none 


Paper: (Clemens et al. 201 1) 


- 


071112C 


7081 


7070 


- 


- 


- 


071117 


7098 


none 


none 


7108 


- 


080303 


7351 


none 


none 


none 


none 


080307 


7362 


none 


7362 


- 


- 


0803 19C 


7442 


7468 


- 


- 


- 


080320 


7473 


none 


none 


7474 


- 


080325 


7531 


none 


none 


none 


none 


0804 13B 


7606 


7598 


- 


- 


- 


080430 


7647 


7647 


- 


- 


- 


080503 


7673 


none 


none 


none 


7673 


080517 


7748 


none 


7748 


- 


- 


080523 


7772 


none 


7772 


- 


- 


0806 13B 


7876 


none 


none 


none 


7873 


080701 A 


7913 


none 


none 


none 


none 


080702B 


7924 


none 


none 


none 


none 


080707 


7952 


7948 


- 


- 


- 


080710 


7969 


7962 


- 


- 


- 


080714 


7979 


none 


none 


none 


none 


080727B 


8030 


none 


8033 


- 


- 


080805 


8059 


8060 


- 


- 


- 


080903 


8176 


none 


none 


none 


none 


080915B 


8234 


none 


none 


none 


none 


0809 16A 


8243 


8254 


- 


- 


- 


081104 


8479 


none 


none 


none 


none 


081121 


8537 


8542 


- 


- 


- 


090129 


8861 


none 


none 


none 


none 


09040 IB 


9068 


none 


none 


none 


9066 


090518 


9393 


none 


none 


none 


none 


090520 


9417 


none 


none 


none 


none 


090530 


9443 


none 


none 


none 


9443 


090726 


9706 


9712 


- 


- 


- 


090727 


9724 


none 


none 


none 


9724 


090904A 


9888 


none 


none 


none 


none 


091018 


10034 


10038 


- 


- 


- 


091024 


10072 


10065 


- 


- 


- 


091026 


10081 


none 


none 


none 


none 


09 1208 A 


10253 


none 


none 


none 


none 


091208B 


10265 


10263 


- 


- 


- 


1001 11A 


10317 


none 


none 


none 


none 


1001 15A 


10325 


none 


none 


none 


none 


100213B 


10412 


10422 


- 


- 


- 


1003 16A 


10501 


none 


none 


none 


none 


100316B 


10500 


10495 


- 


- 


- 


100423 A 


10651 


none 


none 


none 


10658 


1005 14A 


10761 


none 


none 


none 


none 


100702 A 


10926 


none 


none 


none 


none 


100704A 


10929 


none 


none 


10940 


- 


100727 A 


11001 


none 


none 


none 


none 


100802 A 


11031 


none 


none 


none 


none 


1008 14A 


11094 


11089 


- 


- 


- 


100823A 


11135 


none 


none 


none 


none 


1009 17A 


11289 


none 


none 


none 


none 


101008 A 


11318 


none 


none 


none 


11318 


101011A 


11332 


none 


none 


none 


11331 


101017A 


11345 


none 


none 


none 


11345 


101023A 


11363 


none 


none 


none 


none 


101114A 


11405 


none 


none 


none 


none 


101213A 


11448 


11457 








110102A 


11509 


none 


none 


none 


none 
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J.C. Tello et al.: Searching for galactic sources in the Swift GRB catalog 



GRB# 


F.R.E.D. 


Redshift 


"High-z" 


Any z 


Peaks 


11 0223 A 


11764 


none 


none 


none 


none 


Total: 


111 


77 


71 


59 


49 



Table 4. GRBs of the samples used for this study and the GCN used to discrimi- 
nate each burst. 



